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ABSTRACT 

We have obtained medium-low resolution spectroscopy and BVI CCD imaging of Berke- 
ley 32, an old open cluster which lies in the anticentre direction. From the radial velocities of 
48 stars in the cluster direction we found that 3 1 of them, in crucial evolutionary phases, are 
probable cluster members, with an average radial velocity of H-106.7 (cr = 8.5) km s^^. 

From isochrone fitting to the colour magnitude diagrams of Berkeley 32 we have obtained 
an age of 6.3 Gyr, (to - A/)o= 12.48 and E{B - V) = 0.10. The best fit is obtained with 
Z=0.008. A consistent distance, (m — A/)o — 12.6 ± 0.1, has been derived from the mean 
magnitude of red clump stars with confirmed membership; we may assume {m — M)q ~ 
12.55 ±0.1 

The colour magnitude diagram of the nearby field observed to check for field stars con- 
tamination looks intriguingly similar to that of the Canis Major overdensity. 

Key words: Open clusters and associations: general - open clusters and associations: indi- 
vidual: Berkeley 32 - Hertzsprung-Russell (HR) diagram 



1 INTRODUCTION 

Open clusters (OCs) are an important population of the disc of our 
Galaxy and may be used to trace its properties (e.g., spiral struc- 
ture, k inematics, age and age distribution, metallicity, etc; see lFriell 
1 19951 for a review). In particular, by studying the old component 
of the OC population we may be able to trace the disc properties 
at all epochs since its formation, and to derive its evolution. The 
number of known - and well studied - old OCs has steadily grown 
in recent years; what we still miss is a populous sample analyzed in 
a complete and homogeneous way to derive the cluster properties 
like membership, age, distance, reddening, and (detailed) chemical 
abundances using both photometric and low- and high-resolution 
spectroscopic data. 

We have therefore undertaken a project devoted to precise and 
homogeneous derivation of the fundamental properties of a large 
and well chosen sample of old open clusters. The goal of our project 
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is to get insight on the formation and chemical enrichment mecha- 
nisms of the Galactic disc. A description of the motivations of our 
program, and a status r eport of its photometric part are given in 
iBragaglia & Tosi|j200d) . 

As part of this project, we present here a study of Berkeley 32 
(Be 32), an old anticentre open cluster, with coordinates 02000 = 
06''58'"07^ 52000 = +06°25'43", and I = 208°, 6 = -^4.4°. 
We have obtained medium resolution spectra of about 50 stars to 
derive the membership, as we have recently done for Berkeley 29 
jBragaglia. Held & Tosil2005h . while B,V,I CCD images are used 
to derive age, distance, reddening, and approximate metallicity. 

The cluster has already been studied by several authors. 
IScott et all il995h obtained low-resolution spectra and measured 
radial velocities of 14 stars, finding that 10 of them were clus- 
ter members and obtaining an average velocity of 101 km s^^, 
wi th an error of 10 km s^^. These velocities were adopted also 
bv lFriel et all j2002l) . who determined metallicities - with an indi- 
vidual error of about 0.1 dex - and found an average [Fe/H] = — 0.5. 

iKaluznv & Mazuil <199lh obtained U,B,V and Washington 
photometry of Be 32; combining different methods, like the Mor- 
phological Age Ratio (see Sect. 3), the comparison to well stud- 
ied clusters like NGC 188 and M 67, and main sequence fitting, 
they determined that Be 32 has an age of about 6 Gyr, a metal- 
licity [Fe/H] = -0.37 ± 0.05, reddening E{B -V)= 0.16, and 
distance modulus (m - M)v = 12.95 ± 0.15. lRichtler & SagaJ 
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Table 1. Log of the photometric observations. 
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Figure 1. A 20 seconds V exposure on Be 32 (9.4 X 9.4 arcmin-^). The 
stars observed with MOS at DOLORES are indicated by red circles and 
blue boxes for Mask 1 and 2, respectively. 



y001|), working with V, I data, reached similar conclusions; they 
employed isochrone fitting and red clump luminosity to derive an 
age of 6.3 Gyr, [Fe/H] ~ -0.20, E{B - V) = 0. 08, and true dis- 
tance m odulus (m — M)o = 12.60±0.15. Finallv. lHasegawa et alj 
i2004h observed Be 32 (that they indicate with the alternative name 
of Biurakan 8) as part of a survey of 14 anticentre clusters; their 
B,V,I photometry is of somewhat lower quality, and results they 
derive from isochrone fitting and morphological parameter indica- 
tors are not in perfect agreement, but they confirm that Be 32 is 
among the oldest OCs. 

The paper is organized as follows: observations and reductions 
are described in Section 2, the CMDs are presented in Section 3, 
together with results on isochrone fitting; the radial velocities and 
memberships are derived in Section 4; a summary can be found in 
Section 5. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Photometry 

Our data were acquired at the Italian Telescopio Nazionale Galileo, 
on the Canary Islands, using DOLORES (Device Optimized for 
the LOw RESolution), a focal reducer capable of imaging and low 
resolution spectroscopy, on UT 2000 November 26, and UT 2004 
February 14. DOLORES mounts a 2k Loral thinned and back- 
illuminated CCD, with scale of 0.275 arcsec pix~^, and a field of 
view of 9.4 x 9.4 arcmin^. We observed two fields, one centered 
on the cluster (see Fig.0 and a second one about 20 arcmin away, 
to determine the level of field star contamination. 

For each field, TableQlists the date of the observation together 
with the filters used and other relevant information. Seeing condi- 
tions were better for the short exposures taken in 2004; however, 
since this is a rather loose field, even modest seeing did not rep- 
resent a problem for our photometry. We did not have any shutter 
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26 1 1 2000 


04:20:18 
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1.094 


1.93 
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26 1 1 2000 


04:16:15 
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1.091 


2.06 
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26 1 1 2000 


04:14:09 


20 


1.090 


2.48 
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26 1 1 2000 


04:34:15 


300 
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1.93 
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26 1 1 2000 


04:30:20 


120 


1.101 


1.93 


I 


26 1 1 2000 


04:28:00 


20 


1.099 


1.93 



timing problems, since they arise only below 0. 1 s with this instru- 
ment and even our shortest exposure time was 10 times longer. 

The reduction procedure is the same adopte d in the other pa- 
p ers of this series, and deta ils can be found e.g., in lTosi et al.U2004l) 
or lDi Fabrizio et ai]i2005h . Standard IRAF ^ routines were utilized 
for pre-reduction, and the IRAF version of the DAOPHOT-II pack- 
age jstetson 1987, Davis 1994) was used with a quadratically vary- 
ing point spread function (PSF) to derive positions and magnitudes 
for the stars. Output catalogues for each frame were aligned in posi- 
tion and magnitude, and final (instrumental) magnitudes were com- 
puted as weighted averages of the individual values. 

Aperture correction to bring the PSF magnitudes on the same 
scale of the aperture magnitudes, i.e. the same of the standard stars, 
was computed for each relevant frame. The corrections range from 
about —0.15 to —0.35 mag. 

We obser ved the standard areas PG023 1-1-051 and Rubin 149 
jLandolll99l) and derived the calibration equations assuming the 
average extinction coefficients for the site, as given in the site web 
pages (www.tng.iac.es: k_b = 0.25, fty = 0.15, k/ = 0.07): 

B = b + 0.0548 X {b-v) + 1.464 (rms = 0.012) 

V = v- 0.1448 X {b-v) + 1.294 (rms = 0.027) 

V = v- 0.0933 X {v-i) + 1.270 {rms = 0.021) 
I = i + 0.0228 X {v 0.846 {rms = 0.021) 

where b, v, i, are the aperture-corrected instrumental magni- 
tudes, after further correction for extinction and for exposure time, 
and B, V, I are the output magnitudes, calibrated to the Johnson- 
Cousins standard system. 

Finally, we determined our completeness l evel in each ban d 
using extensive artificial stars experiments (see iTosi et aljE004l) . 
whose results are given in Table 2. 

^ IRAF is distributed by the National Optical Astronomical Observatory, 
which are operated by the Association of Universities for Research in As- 
tronomy, under contract with the National Science Foundation 



Table 2. Completeness level for the central and external fields; mag is the 
calibrated magnitude {B, V or /, calibrated with the equations given in the 
text, and assuming b — v and v — i = I). 
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1.00 


1.00 


1.00 


1.00 
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1.00 
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0.85 


0.97 
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0.92 
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0.80 


0.81 


0.78 


0.97 


0.94 


0.87 


18.50 


0.74 


0.69 


0.68 


0.96 


0.93 


0.84 


19.00 


0.67 


0.59 


0.54 


0.93 


0.93 


0.73 


19.50 


0.52 


0.45 


0.37 


0.91 


0.90 


0.52 


20.00 


0.33 


0.28 


0.21 


0.89 


0.86 


0.29 


20.50 


0.16 


0.10 


0.09 


0.85 


0.78 


0.11 


21.00 


0.07 


0.05 


0.03 


0.69 


0.58 


0.04 


21.50 


0.02 


0.01 


0.01 


0.42 


0.32 


0.01 


22.00 


0.00 


0.00 


0.00 


0.22 


0.15 


0.00 


22.50 


0.00 


0.00 


0.00 


0.07 


0.05 


0.00 



We a lso checked out photometry ag ainst the s tudies of 
iKaluznv & Mazur. tl991.) and .Richtler & Sag^ <200ll) . finding 
a very good agreement. The stars in c ommon, about 800 wit h 
iKaluznv & Mazuil Jl99lh and 1700 with lRichtler & SagaJ J200lh . 
respectively, cover the whole interval in magnitude and colours. 
They have a mean difference (our photometry minus literature) 
of AB = -0.003 ((7= 0.070). AV = -0.015 (cr=0.061), for 
lKaluznv&Mazud<199lh and AJ = - ( 3.011 (cr^O.lOl), AV = 
+0.015 (o-=0.072) for lRichtler & Sagail <200ll) (i.e., we are inter- 
mediate between the two photometries in V). No trends with mag- 
nitude or colour are present. 

2.2 Spectroscopy 

The targets for the spectroscopic observations were selected from 
the preliminary photometry of images taken during our first run, 
trying to sample the relevant evolutionary phases: red giant branch 
(RGB), red clump (RC), upper main sequence (MS) and main- 
sequence Tum-Off (TO) point. We also added two stars near 
the RGB tip, taken from published catalogues (stars 2 and 4 in 
IKaluznv & Mazudll99lh . These stars are too bright for the expo- 
sure times of our first photometric images and only one of them 
was later completely recovered in our second run. A few proba- 
ble field stars, based on the preliminary colour-magnitude diagram, 
were also included as a check. The spectra were acquired using 
the Multi Object Spectroscopy (MOS) facility of DOLORES on 
UT February 14, 16, 17 2004; 4 exposures were take with the first 
mask, and 3 with the second one (see Table 3), using the VPH grism 
centered on Ha and a I'.'l slit (resolution 1.36 A, or 7?. ~ 4800). 
The full wavelength range of our spectra is 680 A, but in Fig. |2| 
only a region centred on Hq is s hown. The configurat ion is the 
same already used for Berkeley 29 iBragaglia et all2005h . as is the 
reduction procedure. 

Given the much better S/N of most of the present spectra (see 
Fig.|2|for some representative cases with S/N about 25, 75, 125, 150 
from bottom to top), a different strategy was instead adopted to de- 
rive the radial velocity (RV). To this end we used the task rvidlines 
in IRAF with a list of about 35 lines, mostly of iron, and comprising 
Ha, the only line always visible in all spectra, even at very low S/N 
and/or high temperature. The formal error of the velocities mea- 
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Table 3. Log of the spectroscopic observations 
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Figure 2. A few representative spectra, shifted for clarity; from top to bot- 
tom we show a RC, a RGB and two TO/MS stars 

sured by n'idlines is of about 5 km s~^. The observed velocities 
were corrected to heliocentric, and averaged (see Section 4). 



3 THE COLOUR-MAGNITUDE DL4GRAMS 

The final photometric cluster catalogue^ includes 2031 objects; 
1541 of them are identified in all 3 filters, 1576 at least in B and 
V and 1996 at least in V and /. Similarly, the external field cat- 
alogue has 1398 objects (1031 of them with magnitudes in all 3 
bands, 1069 at least in B and V, and 1348 in V and /). The cor- 
responding CMDs are shown in Fig.|3| where we also anticipate 
results on membership, described in Sect. 4. Pixel coordinates for 
all Be 32 stars were transformed to equatorial coordinates using 
software written by R Montegriffo; we used a linear relation, and 
residuals have a rms of about 0.15 arcsec in both coordinates. 

Some disc star contamination is noticed in the cluster CMD, 

^ The catalog ca n be obtained in electronic form through the BDA 
iMermilUoci 19951) at |http://www.univie.ac.at/webda/new.htinl] 
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Figure 3. (a) and (b) CMD for Be 32 based on our central field. The main sequence and red clump of the cluster are easily discernible. We indicate, only in 
(a) for clarity, the position in the CMD of the stars observed with DOLORES; (red) filled squares and (blue) crosses indicate members and non members, 
respectively, on the basis of their RVs (see Sect. 4). (c) and (d) The CMD for the control field. 



yet we are able to locate the TO and RC at F ~ 16.0, B -V ~ 
0.6, y - / ~ 0.7, and V ~ 13.7, B - V ~ 1.1, 1/ - J ~ 
1.2, respectively. Most stars falling out of the main features are 
contaminating field objects, but there are a few ones, brighter than 
the TO and bluer than the RGB, that may instead be blue stragglers; 
they are worth of further examination. Furthermore, an accurate 
measurement of the RC mean magnitude is obtained using the 4 
RC stars with spectroscopically confirmed membership (see next 



Section). We derive (V> = 13.66 (rm.s. = 0.05) and (J> = 12.53 
(r.m.s. = 0.04), where the errors are the standard deviations of the 
measurements. 

Be 32 appears to be old; a first measure of its age can be ob- 
tained through morp hological age indicato rs, like the Morphologi- 
cal Age Index (MAI: |janes & PhelDsll994l) . based on the difference 
in magnitude between t he TO and the RC {5V), or the Mor pholog- 
ical Age Ratio (MAR: lAnthonv-Twarog & Twaroalll98^ . which 
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Figure 5. The absolute magnitude versus dereddened colour diagrams for: (a) CMa [using (m — M)o=14.5, E{B — V) =0.10 from iBellazzini et all2004l) 1. 
plotting only the fraction of stars implied by the different areas; (b) our external field [using {m — Af)o=13.4, E{B — V) =0.3: this reddening value compares 
very well with the one derived in the Schlegel et al. 1998 maps, which is 0.295 mag in this direction]; (c) the Besan9on model, computed on an equal-area 
field, shifted as our field and adopting the completeness of our data. The lines on the blue of the MS indicate the regions of the CMDs we are comparing. 



also takes into account the difference in colour between the two 
points. From our data we obtain 5V ~ 2.3 and S{B — 1/) ~ 
0.65; these figures have an attached uncertainty of about 0.1 mag, 
since the actual TO and RC points are difficult to locate in a di- 
agram contaminated by field stars (even if we were helped by 
the membership information of selected stars). They translate into 
ages of about 6 Gyr using the IVIAI, and between 5 and 7 Gyr 



using the original jAnthonY -Twarog & Twarog"l983) and the re- 
vised iAnthonv-Twarog & Twarogkl989,.) calibrations of the MAR, 
respectively. Our measure is in good agreement with the age de- 
rived both from morphological age parameters and isochrone fitting 
by Richtler & Sagar (2001). This has to be taken as a first indica- 
tion of the c l uster age, but is substantiated also by the results by 
ISalaris etal]<2004 . They calibrated a similar method, that makes 
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also use of information on tlie metallicity; using the literature val- 
ues of 5V = 2.4 and [Fe/H]=— 0.50, they derived a very similar age 
of 5.91 Gyr. 

We have fitted our data with isochrones, using the models of 
iBertelli et alj il994) . since they come from the same set of evolu- 
tionary tracks that we ha ve already successfully used in our works 
jBraeadia & Tosl2006l and references therein). We tried three dif- 
ferent values of metallicity, Z = 0.02 (i.e., solar abundance), Z = 
0.008 (corresponding to [Fe/H] ~ -0.4) and Z = 0.004 (corre- 
sponding to [Fe/H] ~ —0.7). The first and last sets of isochrones 
did not fit well our data, while isochrones with Z = 0.008 provided 
a good fit. This is in agreement with the m ost recent determ ination 
of the cluster metallicity [Fe/H] = -0.50 <Friel et aljEoO^) .^ The 
results of the best fit, obtained with E{B -V) = 0.10, (m - M)o= 
12.48 and age = 6.3 Gyr (and Z=0.008), are shown in Fig.lH 

An independent estimate of the distance has been derived by 
an empirical comparison of the RC in Be 32 with the RC level 
in open c lusters with similar p roperties and well established dis- 
tance (see lBragaglia et alj2005l for deta ils). We chose as t emplate 
cluster Mel 66, an open cluster for which lSarai edinil i 19991) reports 
[Fe/H] — —0.35 and an age 4.5 Gyr, hence quite similar to Be 32. 
Since the metallicities of Be 32 and Mel 66 are indistinguishable 
within the uncertainties, there is no need to correct for the depen- 
dence of the RC luminosity on metallicity. Instead, we applied a 
small correction for population effects to the V and I luminosity 
of the RC assuming th at Be 32 is 2 G yr older than Mel 66, and us- 
ing the calculations of lGirardi & Salaris (2001.). The correction is 
of the order 0.07 mag (for Z = 0.008) and implies luminosities 
Mv = 0.81 and Mj = -0.27 for the red clump of Be 32. 

With this population correction, and assuming a reddening 
E{B — V) = 0.10 ± 0.05 from our isochrone fit, we obtained an 
extinction-corrected distance modulus (m — M)o= 12.54 ± 0.16 
using Vrc and (m — M)o= 12.62 ± 0.09 from /rc, where the 
errors mainly reflect the uncertainty on the reddening. 

These results do not differ very much from the literature ones 
(see Introduction). We defer any further comparison until a refine- 
ment of our results is obtained using a synthetic CMD technique. 

The CMD of the external field appears to be composed by 
(at least) two populations: a broad young main sequence of field 
stars, crossing diagonally the diagram, and a second more concen- 
trated sequence fainter and bluer than the Be 32 main sequence. 
The latter does not seem, at first sight, completely ascribable to 
the "normal" disc population, as des cribed, e.g., by the Besanjon 
Galactic model bv lRobin et aljj2003l) . A thorough discussion is be- 
yond the purposes of the present paper (and the possibilities of our 
data), but we consider it useful to further check whether we could 
be seeing some extra-disc component. Given the position of our 
field (I = 208°, b — +4.1°) and the distances invo lved we cannot 
have intercepted directly the Monoceros Ring fe.g.. lNewberg et all 
r200Z.Ibata et al. 20ol. This feature of the Galactic disc has also 
been suggested jMartin et aLll2004 to be associated to an over 
density in star counts that could represent the remnant of a dwarf 
galaxy (Canis Major, CMa). While the existence of the Ring itself 
is not disputed, the reality of CMa ha s however been challenged, 
for instance, bv lMomanv et al.l i2004ft who attribute the over den- 
sity to the Galactic Warp. We present in Fig. |5|a comparison of 
our data [Fig. Elb)] to what is thought to be the CMD of CMa 
[ jBellazzini et aL.2Q04) . Fig. (Sja)], and to the modeled disc popu- 



Preliminary analysis of the FLAMES-UVES spectra (see next Section) 
also indicate a metallicity similar to Z = 0.008 



lation [Fig. |5jc)]. All CMDs have been scaled to absolute magni- 
tudes and dereddened colours by correcting for the values given in 
the caption. The apparently better agreement between the CMDs in 
panels (a) and (b) suggest s that we might be see ing a portion of the 
disrupting CMa (see also iBellazzini et al .|2006). 

To say something more definitive, further investigations would 
be required, involving e.g. precise radial velocities of lower MS 
stars, to be compared to the ones that are being derived in the di- 
rection of CMa and of several Ring positions. As indicated in the 
next section, our field stars have an average RV of about 43 km s^^ 
(with an rms of 28 km s~^), similar to that of the disc stars in the 
same direction, as deduced from the Besan^on model (33 km s~^, 
with an rms of 28). However, our velocities refer to stars much 
brighter than those of the faint blue CMD feature which might cor- 
respond to extra-disc components. The bright stars of our sample 
are most likely field objects absolutely compatible with the model 
in Fig. Oc). 



4 RADIAL VELOCITIES: MEMBERSHIP 

We obtained RVs for 48 of the 49 objects observed; one star near 
the TO had very low S/N spectra and we did not determine its RV. 
We averaged the 4 or 3 values for the stars in Maskl and Mask2, 
respectively, after applying the heliocentric corrections. These av- 
erage values are given in Table|3| together with their r.m.s. The RV 
for star 1502 was determined independently for both masks (helio- 
centric RV= 66.92 and 78.37 km s^^ for Maskl and Mask2, respec- 
tively), and we give here the average of the two values. Given the 
precision expected on the basis of the resolving power, the internal 
errors of the RV measures, the standard deviations of the means, 
and this last result, we may safely attach an average error of ± 5 
km s"'^ to our velocities. 

Results of our RV study are summarized in Fig.|6| Panel (a) 
shows all the measured RVs, and the derived mean cluster velocity, 
after elimination of all probable field objects (the limit at RV = 85 
km s~^ is equivalent to a 2.5a clipping). Panel (b) is a histogram 
indicating the mean RV and the tail of field stars. The average clus- 
ter velocity derived from our data is (RV) = 106.7 km s~^ (rms 8.5 
km s^^, 31 stars). 

We have six stars in common with IScott et all il995h . who 
derived RVs with a precision of about 10 km s~^. The RVs are in 
good agreement (see T able |3l, with an average difference (MOS 
minus IScott et alJll995h of about 6 km s ^, that becomes about 3 
km s"'^ excluding the most discrepant star. 

We have used the membership information based on our RVs 
to select RGB and RC stars to be observed with FLAMES-UVES 
(R ~ 40000) at the ESO Very Large Telescope in January 2005. The 
preliminary UVES measurements of the RVs of six stars in com- 
mon are in very good agreement with the present values (P. Sestito, 
private communication): ARV(MOS-UVES) = 2.53 (a = 4.92) km 
s~^, with the difference being dominated by 1 star (without it, we 
have ARV = 0.80, a = 2.78 km s"^). 



5 SUMMARY 

We have presented a photometric and spectroscopic study of stars 
in the old open cluster Berkeley 32. The main results are the fol- 
lowing: 

(i) Analysis of the cluster CMDs confirms that this is an old. 
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Table 4. Radial velocities of stars in Berkeley 32. The first column shows the ID number in our photometry, 
KM and RS indicate the corresponding identifier in Kaluzny & Mazur (1991) and Richtler & Sagar (2001) 
respectively; RA and Dec are the equatorial coordinates at J2000; RV is the hehocentiic radial velocity and cr 
is the standard deviation from the mean RV; RVspj is the radial velocity in Scott et al. (1995). The magnitudes 
come from our catalogue (except for star 2689 - which we measured only in B and / - where B, V come from 
iKaluznv & Mazuj|l99lh . The last column indicates the membership status (C for cluster member, F for field 
stars). The ordering reflects the position of the apertures in each of the two masks. Star 1502 was observed in 
both masks, and the RV given here is the average value. 
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Figure 6. (a) Heliocentric radial velocities for all spectra: filled circles in- 
dicate cluster members (RV > 85 km s^^, and open circles are field stars, 
(b) Histogram of the RVs; red and light blue shadings of different slanting 
indicate the cluster and field stars distributions, respectively. 



: 6 Gyr, 



A/)o= 12.48, 



nearby, rather metal-poor cluster (age : 
V)= 0.10, Z = 0.008). 

(ii) An analysis of the mean luminosity of RC stars provides an 
independent confirmation of this distance to Be 32, yielding (m — 
M)o ~ 12.6 ±0.1. 

(iii) The CMD of the comparison field looks similar to the one 
of CMa, i.e. of the proposed originator of the Anticentre Ring. 

(iv) Membership of 48 stars in the cluster directions has been 
determined on the basis of their RV. 

(v) The cluster average RV has been measured: 106.7 km s^^, 
rms of 8.5 km s^^, based on 31 member stars. 

Further work on this interesting cluster is foreseen, starting 
from the present results: for instance, the information on member- 
ship will be useful when analyzing Be 32 with the synthetic CMD 
technique and the B, V photometry in deriving atmospheric param- 
eters of stars observed with high resolution spectroscopy. 
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